Pragmatic Coarse-Graining of Proteins: Models and Applications.

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

Bacterial 2'-Deoxyguanosine Riboswitch Classes as Potential Targets for Antibiotics: A Structure and Dynamics Study.

The spread of antibiotic-resistant bacteria represents a substantial health threat. Current antibiotics act on a few metabolic pathways, facilitating resistance. Consequently, novel regulatory inhibition mechanisms are necessary. Riboswitches represent promising targets for antibacterial drugs. Purine riboswitches are interesting, since they play essential roles in the genetic regulation of bacterial metabolism. Among these, class I (2'-dG-I) and class II (2'-dG-II) are two different 2'-deoxyguanosine (2'-dG) riboswitches involved in the control of deoxyguanosine metabolism. However, high affinity for nucleosides involves local or distal modifications around the ligand-binding pocket, depending on the class. Therefore, it is crucial to understand these riboswitches' recognition mechanisms as antibiotic targets. In this work, we used a combination of computational biophysics approaches to investigate the structure, dynamics, and energy landscape of both 2'-dG classes bound to the nucleoside ligands, 2'-deoxyguanosine, and riboguanosine. Our results suggest that the stability and increased interactions in the three-way junction of 2'-dG riboswitches were associated with a higher nucleoside ligand affinity. Also, structural changes in the 2'-dG-II aptamers enable enhanced intramolecular communication. Overall, the 2'-dG-II riboswitch might be a promising drug design target due to its ability to recognize both cognate and noncognate ligands.

Shared Binding Mode of Perrottetinene and Tetrahydrocannabinol Diastereomers inside the CB1 Receptor May Incentivize Novel Medicinal Drug Design: Findings from an in Silico Assay.

In recent years, therapeutic compounds derived from phytocannabinoids have brought renewed attention to the benefits they offer to ameliorate chronic disease symptoms. Among cannabinoids, tetrahydrocannabinol (THC) is a well-known component of the Cannabis plant, whose active principles have been studied through the years. Another psychoactive phytocannabinoid, derived from liverworts Radula, perrottetinene (PET), has created interest, especially as a pharmaceutical product and for its legal recreational use. Unfortunately, so far, the interaction mode of these compounds at the type 1 cannabinoid receptors (CB1R) binding site remains unknown, and no experimental three-dimensional structure in complex with THC or PET is available in the Protein Data Bank. Today, many computational methodologies can assist in this crusade and help unveil how these molecules bind, based on the already known pose of a structurally similar compound. In this work, we aim to elucidate the binding mode of THC and PET molecules in both cis and trans conformers, using a combination of several computational methodologies, including molecular docking, molecular dynamics, free energy calculations, and protein-energy network studies. We found that THC and PET interact similarly with the CB1R, in a different conformation depending on the considered diastereomer. We have observed that cis ligands adopted a half-chair conformation of the cycle ring containing the dimethyl group, assuming an axial or equatorial conformation producing a different induced fitting of the surrounding residues compared with trans ligands, with higher interaction energy than the trans conformer. For PET, we have seen that Trp-279 and Trp-356 have a marked influence on the binding. After binding, Trp-279 accommodates its side chain to better interact with the PET's terminal phenyl group, disturbing CB1R residues communication. The interaction with Trp-356 might impair the activation of CB1R and can influence the binding of PET as a partial agonist. Understanding the PET association with CB1R from a molecular perspective can offer a glimpse of preventing potential toxicological or recreational effects since it is an attractive lead for drug development with fewer side effects than trans-THC.

Structural insights into NS5B protein of novel equine hepaciviruses and pegiviruses complexed with polymerase inhibitors.

Infections produced by hepaciviruses have been associated with liver disease in horses. Currently, at least three viruses belonging to the Flaviviridae family are capable of producing a chronic infection in equines: non-primate hepacivirus (NPHV), Theiler's disease-associated virus (TDAV), and equine pegivirus (EPgV). The RNA-dependent RNA polymerases of viruses (RdRp) (NS5 protein), from the flavivirus family, use de novo RNA synthesis to initiate synthesis. The two antiviral drugs currently used to treat hepatitis C (HCV), sofosbuvir and dasabuvir, act on the viral NS5B polymerase as nucleoside and non-nucleoside inhibitors, respectively. Both drugs have shown significant clinical inhibition of viral response. In this work, we aimed to model the NS5B polymerase of the equine hepacivirus (EHCV) subtypes 1 and 2, TDAV and EPgV, to assess whether current direct-acting antiviral drugs against HCV interact with these proteins. Crystal structures of HCV-NS5B were used as templates for modeling target sequences in both conformations (open and closed). Also, molecular docking of sofosbuvir and dasabuvir were performed to predict their possible binding modes at the modeled NS5B polymerase binding sites. We observed that the NS5B models of the EHCV and EPgV shared well-conserved 3D structures to HCV-NS5B and other RdRps, suggesting functional conservation. Interactions of EHCV subtypes 1, 2 and TDAV polymerases with sofosbuvir showed a similar molecular interaction pattern compared to HCV-NS5B, while interactions with dasabuvir were less conserved. In silico studies of molecular interactions between these modeled structures and sofosbuvir suggest that this compound could be efficient in combating equine pathogens, thus contributing to animal welfare.

Studying effects of different protonation states of His11 and His102 in ribose-5-phosphate isomerase of Trypanosoma cruzi: an example of cooperative behavior.

The Trypanosoma cruzi ribose-5-phosphate isomerase B (TcRpiB) is a crucial piece in the pentose phosphate pathway and thus is a potential drug target for treatment of Chagas' disease. TcRpiB residues, such as Cys69, Asp45, Glu149 and Pro47, have confirmed their roles in substrate recognition, catalytic reaction and binding site conformation. However, the joint performance of His11 and His102, in the D-ribose-5-phosphate (R5P) in the catalysis is not well understood. In this work, we probed the influence of different protonation states of His11 and His102 on the behavior of the ligand R5P using molecular dynamics simulations, network analysis and thermodynamic integration. Simulations revealed that a protonated His11 combined with a neutral His102 (His11+‒His102) was able to stabilize the ligand R5P in the binding site. Moreover, calculated relative free energy differences showed that when protonated His11 was coupled to a neutral His102 an exergonic process takes place. On the other hand, neutral His11 combined with a protonated His102 (His11‒His102+), sampled conformations that resembled the catalyzed product D-ribulose-5-phosphate (Ru5P). Network analysis also demonstrated some peculiarities for these systems with some negatively correlated nodes in the binding site for His11‒His102+, and exclusive suboptimal paths for His11+‒His102. Therefore, the combined approach presented in this paper proposes two suitable protonation states for the TcRpiB catalytic mechanism, where an extra proton in either histidines might favor R5P binding or influence isomerization reaction to Ru5P. Our results may guide further in silico drug discovery studies. Communicated by Ramaswamy H. Sarma.

Integrating Molecular Docking and Molecular Dynamics Simulations.

Computational methods, applied at the early stages of the drug design process, use current technology to provide valuable insights into the understanding of chemical systems in a virtual manner, complementing experimental analysis. Molecular docking is an in silico method employed to foresee binding modes of small compounds or macromolecules in contact with a receptor and to predict their molecular interactions. Moreover, the methodology opens up the possibility of ranking these compounds according to a hierarchy determined using particular scoring functions. Docking protocols assign many approximations, and most of them lack receptor flexibility. Therefore, the reliability of the resulting protein-ligand complexes is uncertain. The association with the costly but more accurate MD techniques provides significant complementary with docking. MD simulations can be used before docking since a series of "new" and broader protein conformations can be extracted from the processing of the resulting trajectory and employed as targets for docking. They also can be utilized a posteriori to optimize the structures of the final complexes from docking, calculate more detailed interaction energies, and provide information about the ligand binding mechanism. Here, we focus on protocols that offer the docking-MD combination as a logical approach to improving the drug discovery process.

Thiophenacetamide as a potential modulator to NF-κB: structure and dynamics study using in silico and molecular biology assays.

Nuclear factor kappa B (NF-κB) plays critical roles in the regulation of many pathophysiological processes, including inflammation and immune responses, cell growth and apoptosis. This DNA-binding protein receptor is considered an important molecular target to treat many diseases through host-directed therapy. In this line, several drugs containing thiophene cores have been extensively evaluated due to their ability to interfere on NF-κB translocation to the nucleus. In this work, assays using drug affinity responsive target stability (DARTS) revealed that the parent compound N-(Aryl)-2-thiophen-2-ylacetamide referred to as thiophenacetamide (TAA) specifically binds to the p65 subunit of the NF-κB. Since no experimental binding mode of TAA with p65 is available, we explored TAA within putative sites in silico to gain insights into its possible binding mode and behavior. The binding mode of TAA found in Site 1 formed hydrogen bonds with Lys37 and Asp125 on p65, important residues near DNA-binding region. Molecular dynamics simulations showed the stability of this mode of binding in contrast to the other also tested modes. Our results suggest that TAA binding could occur in regions close to residues responsible for DNA binding, increasing NF-κB protein rigidity and affecting the association between DNA and NF-κB. Communicated by Ramaswamy H. Sarma.